Imprinting Light Phase on Matter Wave Gratings in Superradiance Scattering 
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Superradiance scattering from a Bose-Einstein condensate is studied with a two-frequency pump- 
ing beam. We demonstrate the possibihty of fully tuning the backward mode population as a 
function of the locked initial relative phase between the two frequency components of the pump- 
ing beam. This result comes from an imprinting of this initial relative phase on two matter wave 
gratings, formed by the forward mode or backward mode condensate plus the condensate at rest, so 
that cooperative scattering is affected. A numerical simulation using a semiclassical model agrees 
with our observations. 
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Superradiance from a Bose-Einstein Condensation 
(BEC) offers the possibility of studying a novel physics 
associated with cooperative scattering of light in ultra- 
cold atomic systems. A series of experiments [1-4] and 
related theories [5|-|3 have sparked new interests in mat- 
ter wave amplification [ll, ^ , holographic storage [4] , scat- 
tering spectroscopy in optical lattice coherent imag- 
ing [9], coherent atomic recoil lasing I1Q|412| , and quan- 
tum states storage and retrieval ^IsfOT 

In a typical BEC superradiance experiment, an elon- 
gated condensate is illuminated by an off-resonant pump- 
ing laser pulse along its short axis. Due to the phase- 
matching condition and mode competition, highly direc- 
tional light is emitted along the long axis of the conden- 
sate, in the so-called end- fire modes. Consequently, the 
recoiled atoms acquire a well-defined momentum at ±45° 
angles with respect to the pumping laser direction. These 
atomic modes are referred to as forward modes. This for- 
ward scattering is interpreted as optical diffraction from 
a matter wave grating [1, 5]. Meanwhile, atoms in the 
condensate may scatter photons in the end-fire modes 
back into the pumping mode and recoil at ±135° angles, 
forming the so-called backward modes [3, 7], when the 
pumping pulse is short and intense. This pattern was 
interpreted as a result of diffraction of atoms off a light 
grating. A four wave mixing interpretation was proposed, 
involving two optical fields-the pumping laser field and 
an end-fire mode, and two matter wave modes-the con- 
densate and a mode of momentum [7]. 

There is an energy mismatch of four times the re- 
coil frequency for this backward scattering, due to the 
increased kinetic energy of recoiled atoms [sl, Q. Re- 
cently, a two-frequency-pumping scheme has been imple- 
mented [TBI-ITtI, where the pumping beam consisted of 
two frequency components and the frequency difference 
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FIG. 1. (Color online) Schematic diagram of our experiment. 
The two-frequency pumping beam is incident along the short 
X direction, with a linear polarization along the y direction. 
The atomic side modes are denoted in momentum space, each 
labeled with a pair of integers which describe the order in 
the X and z directions, respectively. Within this notation, 
atoms of the condensate at rest are in mode (0, 0). A forward 
scattering event transfers an atom from mode (n, m) to mode 
(n + 1, m =b 1), and a backward event transfers one to mode 
(n — 1, m =b l).The end fire mode e± is along the long axis of 
the condensate in z direction. 



was controlled to compensate for the energy mismatch 
and excite the backward scattering on a long time scale 
with a weak pump intensity. The presence of the back- 
ward mode in the spectroscopic response [l6| and the 
enhancement of the diagonal sequential scattering [l^ 
have been reported. In those experiments the relative 
phase between the two pumping frequency components 
is maintained constant. Although phase is a very impor- 
tant factor for understanding interference and coherence 
phenomena, phase effects on matter wave gratings and 
cooperative scattering have not been reported so far. 
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FIG. 2. Two- frequency superradiance scattering experimen- 
tal setup. The pump laser beam is split into two beams with 
orthogonal polarizations using a PBS. With two AOMs, each 
beam frequency is shifted about 90 MHz, with a little differ- 
ence of Acj. In order to keep the coherence of the two beams 
after frequency shifting, the phase-locked-loops of the radio 
frequency (rf) sources for the two AOMs are locked to the 
same crystal oscillator. A phase detection circuit is used to 
detect the initial relative phase of beating signal of the two 
frequency beams after they are combined. Note that the beat 
signal of the pump pulse is opposite in phase with the signal 
measured by the photodiode. 



We present here our new experimental results that 
show a possibility to fully control the backward mode 
population as a function of the locked initial relative 
phase between the two pumping frequencies. A theoreti- 
cal analysis of the cooperative process is used to confirm 
and explain this phase dependence. It requires to con- 
sider the relationship between matter wave gratings and 
optical waves beating for the first order scattering. The 
initial relative phase of the pumping beam is imprinted 
into two matter wave gratings, one formed by the con- 
densate with a forward mode, the other formed by the 
condensate with a backward mode. 

In our experiment, a nearly pure Bose-Einstein conden- 
sate of about 2 X 10^ ^'^Rb atoms in the \F = 2, mr = 2) 
hyperfine ground state is generated in a quadrupole-Ioffe- 
configuration magnetic trap, with Thomas- Fermi radii of 
50/im and 5/im along the axial and radial directions [l7j . 
The pumping pulse with two-frequency components is in- 
cident along the short axis of the condensate, with its 
polarization perpendicular to the long axis, as shown 
in Fig. [H This arrangement of polarization induces 
Rayleigh superradiance where all side modes possess the 
same atomic internal state [1]. To get the two- frequency 
pumping beam, a laser beam from an external cavity 
diode laser is split into two equal-intensity beams. Their 
frequencies are shifted individually by acousto-optical 
modulators (AOMl and A0M2) which are driven by 
phase- locked radio frequency signals. Therefore the fre- 
quency difference Au between the two beams can be con- 
trolled precisely, as shown in Fig. [2j After that, they are 
recombined to form our linear-polarized two-frequency 



pump beam, which is red detuned by 27r x 1.5GHz from 
the \F = 2,mF = 2) to 3,m'^ 3) transition. 

The frequency difference Auj is chosen to be 27r x 15kHz 
(the corresponding period T is then 66. 67/is), i.e. four 
times the recoil frequency (4cj^ = 2hkf /M) so that we 
reach the two-photon resonance condition taking into ac- 
count the recoil energies deficit for the backward scat- 
tering. The magnetic trap is shut off immediately after 
the pumping pulse, and the distribution of atomic side 
modes is measured by absorption imaging after 21ms of 
ballistic expansion, as shown in Fig. [3l 

The experiment is repeated with different initial rel- 
ative phases between the two frequency components. 
Monitoring the beating signal cos(Aco't + (/)) on a pho- 
todiode, a light pulse of duration T with a definite initial 
relative phase (j)o = Auto + ^ can be generated. Different 
initial relative phases 0o varying between and 27r, are 
obtained by switching on the pulse at different time to, 
which we call "generation time", as shown in Fig. 3(e). 
The generation time to and the switch off time to -\- T 
can be controlled precisely with an AOM. Figure 3(a) 
and (b) demonstrate the initial relative phases of the 
two frequency pumping beam being 7r/2 and 1.47r, the 
corresponding generation time to being 50/is and 20/is, 
respectively. This phase can be read out related to the 
reference signal corresponding to (/)o = 0. 

In Fig. [3] (c), the backward mode is obvious for a rel- 
ative initial phase (j)o = 7r/2 while in Fig. [3] (d), when 
the phase is <po — 1.47r, backward scattering is almost 
suppressed. The ratio between the backward population 
and the total atom number is plotted in Fig. Inversus dif- 
ferent phases ^o, with a pulse duration equal to 66.67/is 
or one cycle, a time step between two points given by 
Ato = 10/is and each point being the average on four 
experimental data. The error bar is their standard de- 
viation. For each experiment, the atom number in the 
BEG, the initial quantum noise and the temperature are 
not exactly the same. These fluctuations are the main 
reason for the experimental uncertainty. Nevertheless, 
the number of backward scattered atoms shows a high 
sensitivity to the initial relative phase between the two 
frequency components and it is possible to obtain an al- 
most complete cancellation of the backward scattering 
for (j)o = 37r/2. 

Such a result was not expected considering in partic- 
ular the interpretations given in previous experiments. 
Especially, the exponential growth model of scattered 
atoms used in previous articles [1] would not accurately 
account for the phase dependence. Even considering 
the depletion of the original mode of the condensate 
it writes dNs{t)/dt = G{t) {No - Ns{t)) Nsit), where 
Ns is the number of scattered atoms and No the ini- 
tial number of atoms in the condensate [1]. The gain 
G{t) is proportional to the pump intensity I{t). This lo- 
gistic differential equation can be solved analytically as 
Ns{t) = [Noe^^foGit)dty^js^^_i^^NoJ^Git)dt^^ Within 

this model, the number of scattered atoms is determined 
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FIG. 4. (Color online) Experimental results and theoretical 
analysis about the ratio of backward scattered atoms to the 
total atom number vs the initial phase (j)o between the two 
pump beams. Each point is an average on four experimental 
data. The dashed line is a numerical simulation using a semi- 
classical theory depicted in the text with one initial seed and 
c oupling factor ^ = 1.5 x 10^. The dotted line is drawn using 
Eq.® with Aifo,! = 7r/2. 




FIG. 3. Superradiance with different initial relative phases 
and time control sequence, (a) and (b): Initial relative phase 
of 7r/2 and IAtt corresponding to a generation time of 50/is 
and 20/is respectively, (c) and (d): Superradiance patterns 
corresponding to (a) and (b), respectively, (e) Control se- 
quence for generating the pump pulse. Note that there is a 
phase difference of tt between the phase shown in figures (a) 
and (b) and that of the real pumping pulse. 



by the integration of the pump intensity /(t), thus pro- 
portional to /q I{t)dt = Iq Jq [1 + cos{Aujt + (l)o)]dt > 0. 
When the pumping pulse duration equals to an exact cy- 
cle, the integration of the pump intensity is always the 
same, regardless of the relative initial phase between the 
two frequency components. This means, within this non- 
Hnear growth model, the number of scattered atoms is 
independent of the initial relative phase. 

The above model is not valid since more than one scat- 
tered mode should be taken into account. Actually, the 
different modes form matter wave gratings and affect the 
scattering behavior. To give a physical explanation, we 
need to clarify the relationship between the population 
of the different modes and the initial relative phase 00 . 
Here we restrain our analysis to a set of diagonal modes 
: (—1,-1), (0,0), and (1,1). Results are similar for 
the other diagonal modes ((1,-1), (0,0), and (—1,1)) 
because of symmetry. Using the slowly varying enve- 
lope approximation (SVEA) [6], the population of mode 
(—1,-1) is given by its wave function ?/^_i^_i(^, r) : 



1,-1 



(r) 



dt 



J —oo \ 



dt 



1,-1 



dt 



(1) 

We use the semiclassical Maxwell- Schrodinger equations 
to describe the coupled dynamics of matter- wave and op- 
tical fields. Although this theory can not be used to an- 
alyze the initial quantum spontaneous process, the pulse 
in the experiment is about several tens of fis and far be- 
yond its initial quantum characteristic. So the effect of 
the pump beam phase on the initial seeds can be omitted 
to analyze the experimental data [l6|, [l3| • Consider- 
ing the trapped BEG is tightly constrained at its short 
axis degrees and the Fresnel index number of the optical 
field is around 1, we consider a one-dimensional semiclas- 
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sical model including spatial propagation effects. Then 
the envelope function of the end-fire mode optical fields 
with a pump laser consisting of two components having 
same intensity and polarization but different frequencies 
uji = cki and uji — Auj, is given by 




(2) 



where r = 2uJrt and ^ = kiz. ^.(^o) = ^ (l + e^^^^+^o^) 
is connected with the initial relative phase (j)o and 
th e coupling fac tor between light and atom g = 
^JZl^cR| (2ijj^AV). R is the Rayleigh scattering rate of 
the pump component and L the BEC length. It indicates 
that the end- fire mode field e_ is due to the coherence 
between different modes, such as modes (0,0) and (1, 1), 
(—1,-1) and (0,0), and spatial overlap between these 
modes is needed. There exists frequency difference 
between adjacent modes. 

The coupled evolution equations of atomic side modes 
is given by 



dr 



-m*(0o)e-V^o,oe' 



I2r 



(3) 



This equation describes the atom exchange between 
modes (—1,-1) and (0,0) through the pump laser and 
end- fire mode field. An atom in (—1,-1) mode may ab- 
sorb a laser photon and emit it into an end-fire mode 
e_, and the accompanying recoil drives the atom into 
the (0, 0) mode. The atom may absorb an endfire mode 
photon and deposit it into the laser mode, to form the 
backward scattering mode (—2,-2) which process is very 
weak and omitted. Here we have omitted the disper- 
sion term, spatial translation, photon exchange between 
modes and non-diagonal modes connected with e+ such 
as mode (-2,0) and (0,-2) [6, 17]. 

Inserting Eq. ([2]) and ([3]) into Eq. ([T]), we can get the 
evolution equation of A^_i^_i 



g^'-i,-i(r) 
dr 



UJrC 
A2t 



[1 + cos(2r + (^o)] 



[Co,ie^2" + C_i,o] +C.C., 



(4) 



where the envelope function of each side mode V^rn,m 
is written as V^rn,m = |'0m,m|e~*'^"^'"^ with the phase 
^m,m assumed to be space independent. Here we define 

C^',n(r) = ST^diCrnAi^r) = \CrnAr)\e-'''^-^- with 

the phase difference between the matter waves A(/9^^^ = 
9^m,m — V^n,n + ^0,0 — V^- 1,-1 7 and the interference grating 
Cm,n{i, r) = V^o,oV^-i -1 di'^l)m,m'^l^n- Then we have 
A(/Po,i = 2(/9o,o - ^1,1 - V^-1,-1, and A(/p_i,o = 0. 

In this case the above Eq. (|4]) can be expanded using 
trigonometric formula 



S7V_i,_i(r) _ g^L 



{[|Co,i|2cos(2r- A(/:)o,i) 



dr uJrC 

+ cos((/)o + A(/Po,i) + cos(4r + (/)o - A(/?o,i)] 
+ 2|C_i,o|[l + cos(2r + (/)o)]} (5) 



In a first approximation, we replace |Cm,n| by its time av- 
erage \Cm,n\- Equation (5) can then be integrated from 
to one cycle, and the cosine functions of time vanish, leav- 
ing only the cos(^o + ^^o,i) term. A simpler expression 
for N-i^-i at the end of the pumping pulse (Acjt = 27r or 
r = tt) is obtained: 



A/'_i,_i(7r) = -27r— |Co,i| cos((^o + A(po,i) + Q^, (6) 

The first term of Eq. (|6]) indicates that the population 
of side mode (—1,-1) is connected with the module of 
four waves |Co,i| = |V^o,oV^*i -i| /^"^ <^^'|V^o,oV^t,i|, the rel- 
ative phase A(/9o,i and the optical initial relative phase 
(/)o, the coefficient of the cosine function determining the 
amplitude of the oscillation. The second part a approx- 
imated to be a constant in the above expression should 
actually include other terms due to the residual time de- 
pendent part of |Cm,n| and may also slightly depend on 
(j)o . This dependence may be neglected in a first approx- 
imation as the time-dependent terms in Eq. (5) have 
different time-dependent signatures and may cancel each 
other after integration. The main purpose of Eq. ([6j) is 
to show the dependence of the backward scattering on 
the initial relative phase although it is hard for us to 
prove its posit ivity due to the complexity of a. However, 
numerical calculations demonstrate that this backward 
scattered atom number is always positive. 

To give a physical explanation to this initial relative 
phase dependence and because of the assumption that 
the phase (pm,m of each side mode does not depend on 
space, we can focus our attention on the terms Cm,n since 
the spatial integral from — oo to oo in the definition of 
Cm,n does not change the relative phase of the matter 
waves. The term (7o,i = '0o,o'0-i,-i <^^Vo,o'0i,i indi- 
cates the presence of two sequential diffractions by two 
gratings, which we analyze as follows. With the spatial 
coordinates described in Fig. 1, the incident laser light 
traveling in the x direction is diffracted by the grating 
formed by modes (0,0) and (1,1), which is along the 
diagonal of the x — z plane, and can be expressed as 
dz\ilJo,oi^l i\cos{kx -\- kz -\- (pi^i — (po^o)- The scat- 
tered light (end-fire mode light) propagating along the 
long axis of the condensate in —z direction is then 
diffracted again by the grating formed by modes (0,0) 
and (—1,-1), which is also along the diagonal of the 
x — z plane, and can be described as |'^/^o,o'0-i -il cos{kx-\- 
kz — cp-i^-i -\- (po^o)^ resulting in a scattered light travel- 
ing back into x direction which is highly correlated to the 
backward scattered atom number. The spatial integral in 
the former grating describes the propagating effect. Here 
a phase shift of 7r/2 takes place as the integral can be ap- 
proximated to be — I V^coV^i 1 1 sm{kx + /cz -h (fi^i — ^o,o) = 
l'^o,o'0i,il cos{kx-\-kz-\-(pi^i — (/:?o,o + 7r/2). Based on light 
diffraction theory, when a beam is incident at an angle of 
45° onto a grating, which is the combination of two grat- 
ings, the maximum diffraction light along the incident 
direction occurs when these two gratings are in phase. 
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That is, -1 + (/:^o,o = ^1,1 - ^0,0 + 7r/2, or 

TT 

That means A(/po,i = 7r/2. According to our simulation 
we find that the relative phase of the matter waves Acpo^i 
is nearly 7r/2, and there is a little deviation for different 
initial relative phases ^0 of pump beam. The dependence 
of A(/9o,i on refiects the intrinsic nature of nonlinearity 
in the superradiance process since it should not depend 
on the initial relative phase between the two frequency 
components when the pumping pulse duration equals to 
an exact cycle for a linear process. The initial relative 
phase of the pumping beams slightly alters the phase 
matching condition, which maybe results in the diver- 
gence between our simple model and the numerical re- 
sult. However, in our discussed model, the result giving 
A^o,i = 7r/2 is reasonable and meaningful for under- 
standing the scattering picture. 

According to Eq. (j6j) and ([7]) and neglecting a, the ratio 
of backward scattered atoms to the total atom number 
versus the initial phase (j)o is simply a sinusoidal func- 
tion and is plotted as the dotted line in Fig. SI with its 
amplitude adjusted to the experimental data. Backward 
scattering is enhanced when the relative phase is nearly 
(j)o = 7r/2, while suppressed for ^0 close to 37r/2. This 
simple model agrees well with our experimental data ex- 
cept for a small shift of the minimum. 

Furthermore, the ratio of the backward scattering 
atom number to the total atom number is simulated 
numerically using coupled Maxwell- Schrodinger equa- 
tions [il, [13 with an initial seed of one atom in modes 
(1, 1) and (—1,-1), as shown by the dashed line in Fig.|4] 
with the coupling factor ^ = 1.5 x 10^. This simulation 
is in good qualitative agreement with our data but not 
completely satisfactory. The difference may be due to the 
necessity to take into account higher order modes in the 
scattering process and to go beyond ID approximation. 
The matching condition of matter waves may also not be 
exactly satisfied. Describing this mechanism under full 
quantum methods [l8| and considering wave mixing for 
the nonlinear high-order scattering modes is our future 
work. 

To show the features of the position of the maximum 
and minimum of the scattering ratio vs the relative phase 
^0, we simulate it for different pumping factors g in the 
week coupling regime, as shown in Fig. [5l The backward 
scattering ratio increases with the coupling factor, which 
depends on the system parameters like the pump beam 
intensity and detuning, the total number of atoms or the 
BEC shape. The simulation clearly shows that the posi- 
tion of the maximum and minimum of the scattering very 
weakly depends on the experimental parameters, this is 
a generic feature. However, the value of the maximum 
increases with the coupling factor. 

The superradiant scattering is initiated by quantum 
mechanical noise, i.e., spontaneous Rayleigh scattering 
from individual condensate atoms. Subsequent stimu- 
lated scattering and bosonic enhancement lead to rapid 
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FIG. 5. The simulation for the ratio of backward scattered 
atoms to the total atom number vs the initial phase (j)o for 
different weak coupling factors g: 1.0 x 10^ (dashed line, its 
value is amplified by 10 times to show); 1.5 x 10^ (solid line); 
1.6 X 10^ (dotted line) and 1.7 x 10^(dashed-dotted line). 
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FIG. 6. The ratio of backward scattered atoms to the total 
atom number vs the initial phase (j)o between the two pump 
beams for different seeds with ^ = 1.5 x 10^. 

growth of the side- mode populations. The ratio of back- 
ward scattered atoms to the total atom number vs the 
initial phase for different initial seeds in modes (1,1) and 
(—1,-1) are shown in Fig. [6l We can see that the posi- 
tions of maximum and minimum do not vary with differ- 
ent seeds, although the amplitude of the scattering events 
changes obviously. The smaller the seeds, the lower the 
backward scattering amplitude. This shows that these 
features caused by the initial relative phase are not re- 
lated to quantum fluctuation effects. 

The above discussion is for the condition Aujt = 27r, 
which means identical coupling factor between light and 
atoms for the different initial relative phase. If the con- 
dition Aujt = 27T is violated, the integration of the pump 
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FIG. 7. The ratio of backward scattered atoms to the total 
atom number vs the initial phase 00 for the same frequency 
difference 15KHz but different pump times t: 22.22 /xs (dashed 
line); 44.44 fis (dotted line); 66.67 fis (solid line) and 88.89 fis 
(dashed-dotted line). ^ = 1.5 x 10^. 

intensity is different for different relative phases, so that 
the coupling factor is different. The ratio of backward 
scattered atoms to the total atom number vs the initial 
phase 00 is shown in Fig. 7 for the same frequency dif- 



ference 15KHz and different pump times t. From Fig. [71 
we know that the position of the extrema will change at 
different t, while the value of the maximum will increase 
with t. 

In summary, we have studied superradiant scattering 
from a Bose-Einstein condensate with a two-frequency 
pump beam. We have shown that this phenomenon de- 
pends much on the initial relative phase between the two 
optical components. Two matter gratings, one formed 
by the condensate at rest plus the side mode with posi- 
tive momentum, the other by the condensate at rest plus 
the side mode with negative momentum, scatter endfire 
modes photons, which in the end affects the atom scat- 
tering process. We found that the relative optical phase 
is imprinted on matter wave gratings and can enhance 
or annihilate scattering in backward modes. Adjusting 
this phase provides a powerful tool for controlling su- 
perradiant scattering. It is already very beneficial for 
understanding the interaction between matter waves and 
optical waves in cooperative scattering. 
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